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We report equilibrium geometric structures of Cu02, Cu03, CuC>6, and CuOg 1 clusters obtained 
by an all-electron linear combination of atomic orbitals scheme within the density functional theory 
with generalized gradient approximation to describe the exchange-correlation effects. The vibra- 
tional stability of all clusters is examined on the basis of the vibrational frequencies. A structure 
with C s symmetry is found to be the lowest energy structure for CuC>2 while a Y shaped structure 
with C2v symmetry is the most stable structure for CUO3. For the larger CuC>6 and CuO^ 1 clusters, 
several competitive structures exist with structures containing ozonide units being higher in energy 
than those with O2 units. The infra-red and Raman spectra are calculated for the stable optimal 
geometries. 
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I. INTRODUCTION 

The copper-oxygen complexes have received significant 
attention because of their role in various processes such 
as catalytic oxidation, applications in the bio-inorganic 
chemistry, and corrosion to name a few. Numerous stud- 
ies, both theoretical as well as experimental, have been 
reported to elucidate the chemical bonding between cop- 
per and oxygen in such complexes. Most of these studies 
were confined to small sizes of CuO^ systems such as 
CuO and Cu0 2 molecules Q, S II @. These earlier 
studies have been recentlyextended to larger systems by 
Wang and coworkers PJ@] who reported photoelectron 
spectra (PES) studies of CuOa; (x = 1 — 6) species. They 
found that the Cu/O complexes containing either an odd 
or even number of oxygen atoms had similarities among 
themselves with respect to their PES. They also found 
several new copper-oxide molecules including CuO^ and 
proposed some structures for these molecules on the ba- 
sis of their PES study. Massobrio and coworkers 0, 
subsequently performed a systematic search for the low- 
est energy structure for CuO^- In their density functional 
calculations with ultrasoft pseudopotentials they found 
the lowest energy structure of CuC>6 to possess ozonide 
Cu(C>3)2 units contrary to the prediction in RefQ of an 
oxide structure Cu(C>2)3 on the basis of photo-electron 
spectra. They also found that the structures containing 
the Cu(C>2)3 units to be roughly within 0.2 eV. Cu is 
known to form OCuO or Cu(Oa) upon reaction with O2 
only under UV radiation at low temperature . The Cu 
atom at room temperature can form only Cu(02) com- 
plexes with O2 because the O atom transfer reaction is 
endothermic. 
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The present work is aimed at providing a detailed pic- 
ture of the stationary points on the potential energy sur- 
face of the CuC"2, CuC"3, CuCe and anionic CuOg clus- 
ters. For this purpose, we perform a systematic search 
for possible candidate structures for the ground state ge- 
ometries of these clusters, followed by calculation of the 
vibrational spectra to identify the nature of stationary 
points. As the energy differences between the isomers 
are small, we employ an all-electron treatment thereby 
eliminating any question due to the ultrasoft pseudopo- 
tential approximation. The present work complements 
the previous investigation on these systems by Pouillon 
and coworkers HHlg. Additionally, we have also calcu- 
lated vertical detachment energies of the CuOe isomers 
and also the infrared absorption and Raman scattering 
spectra of the lowest energy CuC^ isomers. These may be 
useful for comparison with experimental measurements. 

In the following section the computational procedures 
employed in this work are briefly described followed by 
the results and discussion in the subsequent section. 



II. DETAILS OF COMPUTATIONAL 
PROCEDURE 

The calculations have been performed within the 
Kohn-Sham formulation of the density functional the- 
ory with the generalized gradient approximation (GGA) 
to describe the exchange-correlation effects . We have 
used the PBE functional by Perdew and coworkers 
Structural optimizations of all the clusters were per- 
formed by starting from a number of possible candidate 
structures, including those already reported in the lit- 
erature 0, 0|- The optimization was performed using 
the conjugate-gradient as well as the LBFGS method. 
For the case of the anionic CuOe systems, the optimized 
neutral geometries are used as the starting geometries. 
All calculations have been performed using the Naval 
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Research Laboratory Molecular Orbital Library (NRL- 
MOL) package Q- This suite of programs employs a 
Gaussian basis set where the exponentials of the Gaus- 
sian basis are optimized for each atom |10| . The integra- 
tions are carried out using a grid based method with a 
variational mesh. Further, the magnetic moment is op- 
timized during the self-consistency cycle from electron 
counting considerations. The geometries were optimized 
until the forces were less than 0.001 a.u.. The nature of 
stationary points (optimized geometries) on the potential 
energy surface was examined by computing the vibra- 
tional frequencies. These were calculated by displacing 
the atoms from their equilibrium positions by ±0.07 a.u. 
along each Cartesian direction and constructing the dy- 
namical matrix from the forces of the distorted molecules 
The vibrational frequencies were subsequently ob- 
tained by diagonalizing the dynamical matrix. This also 
allows one to compute the zero-point energies. We also 
calculate the IR and Raman spectra of the optimized 
CuOe molecule which are obtained from the derivatives 
of the dipole moments and polarizability tensors |2*7j . 



III. RESULTS AND DISCUSSION 

Our test calculations performed for O2 and O3 
molecules yield the bond lengths of O2 and O3 as 1.22 and 
1.28 Awhich are in excellent agreement with the experi- 
mental values of 1.21 and 1.28 A. The bond angle for O3 
is also well predicted, 118° (calc.) vs. 117° (expt. [^h 
The calculated vibrational frequencies for O2 and O3, 
1547 cm _1 and 708, 106L1185 cm _1 agree well with the 
experimental values |28|,|29j of 1580 cm _1 and 705, 1042, 
1110 cm~ 1 respectively. However, the calculated dissocia- 
tion energies of O2 and O3 are 6.80 eV and 9.06 eV while 
the experimental values are 5.1 and 6.2 eV. Although the 
dissociation energies are grossly overestimated, there is 
good agreement for structural parameters between the 
calculated values and experiment. In the following we 
discuss the CuO^ complexes for x=l,2,3,and 6. 



A. CuO 

The calculated values of bond length, dissociation en- 
ergy and the vibrational frequency of the CuO molecule 
are given in Table [I] Earlier theoretical values calculated 
at different levels of theory along with the experimen- 
tal values are also included for comparison. Our calcu- 
lated value of bond length 1.72 Ais in excellent agreement 
with experimental value of 1.73 A[2j|. The PBE func- 
tional seems to provide significant improvement over the 
LSDA and also seem to perform better than the B3LYP 
in the bond length prediction. Our calculated value of 
the vibrational frequency (680 cm -1 ) is again in good 
agreement with experiment (631 cm -1 ). The dissocia- 
tion energy is however overestimated in all the theoretical 
calculations including the present one. 



TABLE I: The calculated values of bond length (R e in A), 
dissociation energy (D e in eV), and the frequency (w in cm -1 ) 
of CuO are compared with the experimental values as well as 
earlier theoretical values. 



Re 




Id 


EA 


Method 


A 


eV 


cm -1 


eV 




1.73 


2.75 


631 


1.78 


EXPT.° 






631 




EXPT. 6 


1.78 




640 




EXPT. C 


1.72 


3.54 


680 


1.84 


PBE d 


1.74 




592 




CASSCF e 


1.73 


3.58 


630 




PWPP-GGA / 


1.68 


4.13 


670 




PWPP-LSDA-^ 


1.77 




586 




B3LYP 9 


1.69 


4.01 


735 


1.68 


LSDA h 




"Experimental results, Ref. 
''Experimental results, Ref. 
'Experimental results, Ref. 

''Present Calculations (Zero-point energy is included in _D 
'Complete active space self-consistent field method, Ref. 
•'Plane wave pseudo-potential method, Ref. 
fRef. 
ft Rcf 



B. C11O2 

The optimized structures of various isomers of the 
Cu02 molecule are shown in Fig. ^ 

The bent CuOO structure (A in Fig^l with Cs sym- 
metry is the lowest energy structure amongst the Cu(02) 
complexes. It is lower than the linear OCuO molecule 
(structure B) by 0.76 eV. The third structure with C21, 
symmetry ( C in Fig. ^ ) in which the Cu atom occu- 
pies a side on position with respect to the O2 is higher 
in energy by 0.55 eV. The fourth isomer with C2„ sym- 
metry (isomer D in Fig. ^) is not reported in earlier 




FIG. 1: The various isomers of the CUO2 cluster. The isomer 
A is the lowest energy structure. Isomers A, B, and D are 
vibrationally stable. 
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TABLE II: Calculated bond lengths (R e in A), atomization energies (D e in eV), and vibrational frequencies (w in cm of 
Cu02- The experimental values as well as earlier theoretical values are included for comparison. 



Structure 


Method 


Re 




ZOCuO /is 


D e 




LU 2 


^3 






Cu-0 


O-O 






cm -1 


cm -1 


cm -1 


Bent (A) 


DFT/PBE" 


1.88 


1.29 


1 


2.57 


240 


452 


1200 




DFT/B3LYP / 
CI d 


1.91 


1.3 






217 


438 


1188 




1.88 


1.34 


120 1 












cr 


1.91 


1.34 


106 1 












CP 






1 




220 


535 


1161 




DFT/B3LYP* 


1.92 


1.3 






214 


439 


1175 




CASSCF 9 










96 


503 


1359 


Linear (B) 


EXPT." 










193 


658 


819,823 




DFT/PBE 6 


1.67 




1 


2.32 


146 


705 


904 




DFT/GGA C 
CI d 


1.68 






2.27 


122 


663 


864 




1.77 




1 












CP 


1.74 




1 




591 








DFT/B3LYP / 

DFT/PBE 6 

DFT/B3LYP / 


1.65 








142 


727 


953 


Side-on (C) 


1.97 


1.39 


41 1 


2.39 


287i 
171 


260 
401 


940 
1036 




DFT/B3LYP* 
CI d 


1.97 


1.37 


70 




173i 


453 


1119 




2.08 


1.49 


42 












CASSCF 9 


1.99 


1.86 






224 


391 


1052 




EXPT/ 












540,548 


1089 


Wide-angled (D) 


DFT/PBE 6 


1.73 


3.39 


157 3 


2.40 


175 


573 


730 



"Experimental results from Ref. ^3 and lld 
Ref. Pr esent Calculations, Zero-point energy is included in D e 
Ref. 1^ 
21 
22 



d Ref. 
e Ref. 
fRei. 
sRef. 
4 Ref. 
'Ref. 
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works and is found to be lower than the structure C 
by 0.02 eV. In this wide-angled isomer the distance be- 
tween the two oxygen atoms is much larger than that in 
the third structure. The OCuO angle in C is 41° while 
that in the isomer D is 157°. Our vibrational frequency 
calculations have shown that the side on C isomer has 
one imaginary frequency. Such imaginary frequency for 
this isomer has also been observed in the earlier density 
functional based work of Barone and Adamofbl. Earlier 
quantum chemical calculation by Bauschlichcr 23J found 
isomer C to be more stable than the bent structure A. 
They however note that the two structures are energeti- 
cally very close with their order being dependent on the 
numerical details. More recent quantum chemical calcu- 
lation (CCSD) using the effective core potential (ECP) 
reports the side-on structure C to be the lowest. The 
differences are possibly due to the use of ECP. The bent 
CuOO structure (A) found as the lowest energy struc- 
ture in the present work is in agreement with the earlier 
DFT calculations However, the order of isomers as 
the structures A, C, and B in increasing energy ob- 
served in earlier DFT works is now altered in the present 
work due to the new isomer (D) found in this work. This 
isomer is the next low-lying isomer to the lowest energy 



structure. The calculated values of the Cu-0 and O- 
O bond lengths, ZOCuO, spin magnetic moment, atom- 
ization energy and vibrational frequencies are compared 
with the existing theoretical and experimental results in 
Table [n] The earlier theoretical calculations have pre- 
dicted the dissociation energy of Cu02 with respect to 
Cu and 2 to lie between 0.42 -0.56 eV GUM Eg. The 
lowest-energy structure has a doublet state in agreement 
with the experimental observation [eH ITfjl ] . Our results 
show that in the linear structure B, the Cu-0 distance 
is 1.67 Awhich is in good agreement with the prediction 
based on the experimental absorption spectra|18j and is 
also consistent with the findings of another experimental 
result 0. Another density functional based calculation 
flU using the Perdew-Wang91 GGA 26], however pre- 
dicts the quartet state to be lower in energy by 0.1 eV. 
Since the two GGA functionals, PBE in the present and 
PW91 in Ref. 0], are essentially similar, the prediction 
in Ref. 0] could be due to other numerical details like 
choice of basis set in Ref. They however conclude the 
ground state to be a doublet on the basis of vibrational 
frequencies. Density functional calculation by Cherithin 
et al. 0] using the hybrid B3LYP GGA also predicts the 
quartet to be lower than the doublet, but again assign the 
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doublet to be the ground state on the basis of frequency 
calculations. Our predicted values of frequencies are in 
good agreement with the experimental values and earlier 
B3LYP calculation^^. Unlike previous theoretical cal- 
culations, the present calculation correctly points doublet 
to be the ground state on the consideration of both the 
dissociation energy as well as the frequencies. 

The frequency labeled u)± in Table II corresponds to a 
bending mode in the bent structure or a flapping mode 
for the linear and the wide-angled structures. The wi 
refers to a stretching of Cu-0 bonds resulting in a breath- 
ing mode of vibration for the symmetric structures. The 
third frequency corresponds to the 0-0 stretch mode in 
the bent structure while in the linear and wide-angled 
structures it corresponds to the vibration of the Cu atom. 
The 0-0 stretching mode of the bent structure at 1200 
cm -1 shows a large peak in the infrared absorption spec- 
trum while the mode at 240 cm is strongly Raman ac- 
tive. 



C. CuO :i 

The equilibrium structures of Cu03 are presented in 
Fig. [3 The first structure (A in Fig. |2J) has a Y shape 
with Cu atom occupying the center of the fork position. 
The isomer B has a distorted rhombus structure while C 
and D structures are chain-like. One of the chain struc- 
tures has a linear O-Cu-0 unit while the other one has a 
wide-angled O-Cu-0 unit. The spin magnetic moments, 
dissociation energies and the vibrational frequencies of 
the four isomers are given in Table IIIII Our results in- 
dicate that the Y shaped cluster containing a side-on 
Cu02 unit is the lowest energy structure. The spin state 
of this structure is quartet. This result is in agreement 
with car licr published DFT calculations HQS- The dis ~ 
torted rhombus structure has the Cu atom bonded with 
an O3 unit which we refer to as copper ozonide structure. 
This isomer is higher in energy by 0.54 eV. The chain- 
like structures have energies lying between the Y shaped 
structure and the ozonide structure. 

The calculated vibrational frequencies of all the four 
isomers are real indicating that all the four structures are 
local minima on the potential energy surface (Cf. Table 
IIII|) . In case of the lowest energy structure (Structure A), 
the lowest two modes correspond to bending motions, the 
first one in plane and the second one perpendicular to the 
plane of the molecule. The third mode pertains to the 
rotational motion of the O2 unit about an axis perpendic- 
ular to the 0-0 bond. The next two modes correspond to 
Cu-0 bond stretch modes. The largest frequency (1158 
cm -1 ) mode refers to the 0-0 stretch mode. Similar to 
the smaller clusters, this mode of vibration is strongly IR 
active while Cu-0 stretch modes show significant Raman 
scattering. 

In the ozonide isomer (B), the lowest three frequencies 
respectively involve the asymmetric stretch of the Cu-0 
bonds, vibration of the central O perpendicular to the 




FIG. 2: The various isomers of the CUO3 cluster. 



plane of the molecule and a symmetric stretch of the Cu- 
O bond leading to a breathing mode. The frequencies 
from 687 to 987 cm -1 pertain to the flapping , asymmet- 
ric 0-0 stretch and a symmetric 0-0 stretch modes of 
the ozone unit. However, the frequencies are lower than 
a free ozone unit indicating weaker 0-0 bonds in the 
O3 unit due to the presence of the Cu atom. The mode 
at 807 cm -1 shows strong infrared absorption while those 
around 200 cm _1 are highly Raman active. Cherithin et 
al. 0| have observed peaks around 800 cm -1 (812, 802 
and 792 cm -1 ) in their IR measurement which they have 
attributed to the Cu-ozonide structure. The present re- 
sults are supportive of this assignment. 

In the two chain-like structures the vibrational fre- 
quencies are comparable and the modes are same. The 
lowest two modes correspond to the vibration of the Cu 
atom in and out of the plane of the molecule. The next 
two frequencies refer to the asymmetric Cu-0 stretch 
modes while the largest frequency corresponds to the O- 
O stretch mode which is also strongly IR active. The vi- 
brational frequency of the 0-0 unit in the linear modes 
tends towards that of the free O2 frequency (1547 cm -1 ). 
The Cu-0 stretch modes display Raman activity. The 
dissociation energies of the chain- like molecules C and D 
into CuO and O2 are 1.14 and 0.76 eV. 
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TABLE III: The spin magnetic moments (in bohr magneton), 
atomization energies (eV) and vibrational frequencies of the 
various isomers of CUO3. 



Isomer Method 



[IB D e OJe Cm 



EXPT. a 

A DFT/PBE C 3 2.92 
DFT/PW91 d 3 2.90 
DFT/B3LYP 6 3 

B DFT/PBE C 1 2.79 
DFT/PW91 d 1 
DFT/B3LYP 6 1 
EXPT. 6 

C DFT/PBE C 3 2.87 
DFT/PW91 d 3 
DFT/B3LYP 6 3 

D DFT/PBE C 1 



660(90) 

136, 143, 273, 378, 703, 1158 

366,683,1160 

386, 656, 1230 

195,229,321,698,807,987 

697, 799, 977 

724, 879, 1055 

792, 802, 812 

69, 95, 206, 458, 676, 1267 

441, 654, 1258 

400, 635, 1351 

98, 163, 277, 468, 800, 1358 



"Experimental results, R.cf. 2. 
6 Ref. 

c Ref. Pr esent Calculations, Zero- 
d Ref. [3 



point energy is included in De- 



li. Cu0 6 

In the C11O3 clusters, we have seen that the ozonide 
cluster is higher in energy than the Y shaped structure 
which contains an O2 unit. On the other hand, the CuC>2 
triangular structure is higher in energy than a linear Cu- 
0-0 structure. It would be interesting to see how the 
cluster structure changes in the presence of more oxygens. 

The knowledge of structures for the smaller sizes is 




FIG. 3: The various isomers of the CuOs cluster. The isomer 
A is the lowest-energy isomer. 



useful in building up the starting geometries for the op- 
timization of CuOg- We have started with several possi- 
ble structures as input for optimization. The optimized 
geometries are presented in Fig. [3J The first two struc- 
tures, labeled as A and B in the FigO contain the three 
O2 units and also a Cu02 triangle. The difference in 
the structures lie in the orientations of the open-ended 
O2 units. The isomers C and D contain closed CU-O3 
rings. The isomer C is planar and the isomer D is 3- 
dimensional. The structure E has both O3 and O2 units 
of which the O3 forms a closed ring with the Cu atom. 
For the next isomer (F ) also, both O3 and O2 units 
are present and the Cu atom forms end-on bonds with 
these units. The last structure G has three O2 units 
all of which form end-on bonds with the Cu atom. This 
structure was arrived at by starting the geometry op- 
timization with the structure predicted by Wang et al. 
based on their photoelectron spectra. Our calculations 
have shown the structure labeled as A has the lowest 
energy among all the structures shown here. The struc- 
tures B and G are comparable in energy and both are 
higher relative to structure A by 0.12 and 0.22 eV re- 
spectively. The isomers C and D with ozonide subunits 
are higher than A by 0.95 and 0.94 eV respectively. The 
ozonide planar structure C was reported to be the lowest 
energy structure by Pouillon and Massobrio Q ■ This dif- 
ference in the present and their calculations for the lowest 
energy structure is probably due to their choice of ultra- 
soft pseudopotential for representation of core electrons 
and brings out the importance of treating all electrons at 
equal footing in the CuO^ system. The structures E and 
F are higher in energy by 1.62 and 1.79 eV, respectively. 
The atomization energies of these isomers are presented 
in Table IIVI The atomization energy includes the zero- 
point energy of the systems. 

The vibrational frequencies of the isomers are pre- 
sented in Table IIVI Both the structures A and B are 
found to be vibrationally stable. In these two lowest 
energy structures A and B , the modes of vibrations 
are pretty similar. The low frequency modes are asso- 
ciated with the bending and flapping of the open-ended 
O2 units. On the other hand the high frequency modes 
are associated with the stretching of the 0-0 bonds re- 
flecting the strength of these bonds. The Cu-0 bond 
stretch modes have frequencies lying between 250 to 525 
cm -1 . The lower frequencies mainly correspond to the 
bond bending motions of the O2 units. In the struc- 
ture A the high frequency modes corresponding to the 
0-0 bond stretching show significant IR absorption while 
the modes below 400 cm -1 show Raman activity. On the 
other hand, the structure B is IR inactive and has modes 
below 200 cm _1 that are Raman active. 

The structure C comprising an ozonide unit is also vi- 
brationally stable. This structure is both Raman and IR 
active. The Raman scattering is seen in the low frequency 
range of 300-350 cm -1 , where the modes correspond to 
bending of the molecule and also symmetric and asym- 
metric stretching of the four Cu-0 bonds. The IR active 
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modes are the high frequency modes in the range of 800 
- 920 cm -1 . These are again 0-0 stretch modes which 
are asymmetric. The largest mode is also a symmetric 
0-0 stretch mode leading to a breathing mode. How- 
ever, unlike the other structure where the largest mode 
is strongly IR active, this mode is IR inactive. This is 
due to the fact that this mode is symmetric and does not 
lead to a change in dipole moment of the molecule. 

The structure D with two orthogonal planar O3 units, 
is a saddle point with an imaginary vibrational frequency. 
The eigenvector of the mode shows it related to a bending 
motion which try to break the orthogonality of the two 
ozone units. This orthogonality was not enforced during 
the optimization. Both the structures E and F contain- 
ing both O2 and O3 units are found to be vibrationally 
stable. However, they are very high on the energy scale. 
This is again consistent with the energy ordering seen in 
the CUO3 isomers in which both the ring-like and linear 
structures are higher in energy than the Y shaped one. 
The vibrational analysis of the structure G showed it to 
be a saddle point on the potential energy surface. The G 
isomer containing three open-ended O2 units shows IR 
activity in the range of 1200-1250 cm -1 . These again are 
0-0 stretching modes. However the largest frequency 
mode which corresponds to symmetric stretching of all 
the 0-0 bonds is found to be IR inactive. This can be 
attributed to the fact that this mode retains the overall 
symmetry of the molecule and does not lead to a change 
of the dipole moment. 




FIG. 4: The isodensity surface of the lowest energy CuC>6 
structure. 



We have also considered various possible anionic struc- 
tures of the CuOg molecule. The optimized structures of 
the anion are very similar to their neutral counterparts 
and are therefore not shown. These structures are ob- 
tained by starting the optimization with the geometries 
of neutral clusters. As in the neutral case, the struc- 
ture A is found to be the lowest-energy structure. The 
relative energies, vibrational frequencies, vertical detach- 
ment energies of these structures are given in Table [S 
An analysis of the vibrational frequencies show that the 
structures A, B, C and D pertain to local minima on 
the potential energy surface. The structure E is the only 
structure to become somewhat modified from its neutral 
geometry in that the oxygen dimer tilts away from the 
Cu and CU-O2 bond length increases to 2.3 A. Also the 
CUO3 closed ring structure becomes opened at one end. 
This structure is vibrationally stable but energetically 
high. On the other hand, the high energy structure G 
has two imaginary frequencies associated with it. The op- 
timization of the anionic counterpart of structure F was 
abandoned due to failure to obtain a convergent solution. 
The structure D in neutral form is vibrationally unstable 
but the anionic form is vibrationally stable. The adi- 
abatic electron affinity of the lowest energy structure is 
found to be 2.75 eV in which the zero point motion of the 
molecule is also taken into consideration. The calculated 
values of the VDE (Cf. Table© are in the range 2.8-3.5 
eV. The experimental value of the electron binding en- 
ergy in CuOe as measured from photoelectron spectra is 
3.2 eV [2|. The structures with ozonide units as well as 
with O2 units, both have the VDE around this value and 
based on the consideration of the VDE alone, it is diffi- 
cult to infer about the correct anion geometry. However, 
the VDE combined with total energy values indicate the 
oxide structure to be the candidate for the ground state 
geometry for the neutral CuOg- 

The analysis of nature of bonding based on the charge 
density differences shows the expected behavior, that the 
0-0 bonding is covalent while Cu-0 bonds are more ionic 
in nature. The plot of isosurface of charge density are 
shown in Fig. 21 for the case CuOq. The lowest energy 
CuOe is stable by 2.72 eV against dissociation into Cu 
atom and three O2 molecules. This translates to about 
0.68 eV per Cu-0 bond which is a relatively weak bond. 
A simple analysis of integrating charge density around 
each atom inside an atomic sphere of arbitrary radius 
shows the Cu atom to be in the +1 charge state and the 
O atoms in the -2 charge state. There is very little spin 
density on the Cu atom and most of the spin density is 
seen on the O2 units. This is consistent with the param- 
agnetic nature of O2 molecule. The spin density is largest 
on the O atoms farthest from the Cu atom. 



IV. SUMMARY 

In conclusion, we have performed a detail examina- 
tion of various possible structures of the CuO^ molecule 
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TABLE IV: The atomization energies, spin states and the vibrational frequencies of the CuC>6 isomers. The dissociation energies 
include zero-point energy. 



Isomer 


D e (eV) 


fJ-B 


uj e cm 1 


A 


3.29 


3 


30, 95, 107, 128, 148, 174, 264, 294, 299, 405, 476, 496, 1147, 1250, 1324 


B 


3.27 


1 


36, 69, 93, 101, 124, 141, 192, 236, 359, 381, 439, 525, 1168, 1233, 1298 


C 


3.15 


1 


70, 94, 116, 311, 312, 347, 365, 437, 456, 706, 722, 810, 826, 916, 1015 


D 


3.16 


1 


64i, 49, 88, 167, 246, 295, 342, 364, 389, 667, 751, 789, 909, 912, 1030 


E 


3.06 


3 


53, 84, 103, 157, 170, 232, 303, 387, 494, 636, 695, 790, 993, 1340 


F 


3.03 


1 


29, 74, 87, 104, 128, 155, 183, 299, 364, 421, 570, 703, 752, 1179, 1400 


G 


3.26 


5 


25i,44,45,78, 94, 104, 175, 184, 221, 373, 387, 398, 1240, 1251, 1323 



TABLE V: The magnetic moment M (in Bohr magneton), relative energies AE, vertical detachment energy (VDE) and the 
vibrational frequencies of the CuOe anionic isomers. 



Isomer MQjb) AE (eV) VDE (eV) co e cm' 1 

A 2 MOO 3T2 51, 71, 124, 130, 158, 185, 248, 265, 398, 413, 491, 499, 1138, 1209 

B 2 0.06 3.35 27,103, 131, 179, 185, 226, 249, 257, 417, 447, 489, 546, 907, 1133, 1201 

C 2 0.53 3.25 58, 101, 127, 261, 286, 315, 324, 413, 458, 647, 669, 729, 749, 831, 939 

D 2 0.49 3.33 29,79,80, 233, 234, 238, 291, 291, 397, 665, 794, 794, 932, 946 

E 2 0.89 3.49 20, 37, 62, 80, 105, 114, 172, 220, 232, 410, 507, 681, 707, 1143, 1299 

G 4 0.36 2.80 63i, 31i, 51, 65, 105, 130, 166, 196, 204, 379, 390, 405, 1093, 1105 



x = 1—3, 6. The vibrationally stability of the isomers was 
examined by carrying out an analysis of the vibrational 
modes. A bent structure with C s symmetry is found to be 
the most stable structure for CuC>2. On the other hand, 
a Y shaped structure with C2 V symmetry is the lowest 
energy structure for CUO3. Our study shows the poten- 
tial energy surface is rather fiat for CuOg- Within our 
fairly exhaustive search we find two energetically nearly 
degenerate isomers as possible candidate structures for 
the ground state. Both of these structures contain Cu02 
subunits. This observation is consistent with the argu- 
ment by Wang and coworkers based on the PES but it 
rules out their proposal of the structural form. Further, 
contrary to the earlier pseudopotcntial based study, the 
structures with ozonidc subunits are found to be high 



energy isomers. The anionic clusters show trends similar 
to those of their neutral counterparts. The CuOg cluster 
in its most stable form is found to be strongly IR ac- 
tive. The IR active modes correspond to the 0-0 bond 
stretching and generally have high frequencies. 
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